Multi-magnon bound states in easy-axis ferromagnetic zigzag spin chain 
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(Dated:) 

The frustrated spin- 1/2 chain with weakly anisotropic ferromagnetic nearest-neighbor and an- 
tiferromagnetic next-nearest-neighbor exchanges is studied. We focus on the excitation spectrum 
and the low-temperature thermodynamics in the ferromagnetic region of the ground state phase 
diagram. It is shown that the excitation spectrum of the model is characterized by the existence 
of the multi-magnon bound states. These excitations determine the low-temperature magnetic sus- 
ceptibility. The energy of the bound magnon complexes is found using the scaling estimates of 
the perturbation theory and numerical calculations. The relation of the considered model to the 
' edge-sharing cuprate Li^CuOa is discussed. 
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I. INTRODUCTION 



The quantum spin chains with nearest-neighbor (NN) J\ and next-nearest-neighbor (NNN) interactions J 2 have 
been a subject of numerous studies [l|. The model with both antiferromagnetic interactions Ji, J 2 > (AF-AF 
model) is well studied @, [3, 0, S @, 0| ■ Lately, there has been considerable interest in the study of F- AF model with 
the ferromagnetic NN and the antiferromagnetic NNN interactions (Ji < 0, J 2 > 0) @, H Hfl, [U QJ, • One of 
the reasons is understanding of intriguing magnetic properties of a novel class of quasi-one-dimensional edge-sharing 
. copper oxides, which are described by the F-AF model jll, ITU. liH fl7l flM . 
The Hamiltonian of the spin- 1/2 F-AF model is 

N N 

H = J X Y,(SZS* +1 + SlSl +l + Ar5^ +1 ) + J 2 £(S*S* +2 + S*S*^ 2 + A 2 S*S* +2 ) (1) 

n— 1 n— 1 

O 

, where Ji < and J 2 > 0. 

The isotropic case of this model (Ai = A 2 = 1) is intensively studied last years [ll|, [IH, [2(| HH, [HJ . The model with 
y—i • the anisotropy of exchange interactions is less studied. Though the anisotropy in real chain is weak [23[, it influences 
on the properties of such compounds. The ground state phase diagram of model ([I]) with small anisotropy has been 
studied by us in pMj . As shown in Fig.l, the phase diagram consists of three phases: commensurate spin- liquid gapless 
' phase, the incommensurate phase with spin correlations of a spiral type and fully polarized ferromagnetic phase. 
Though this phase diagram is related to the particular case of the anisotropic NN and isotropic NNN interactions, 
• ■ model ([T]) with both A1 ^ 1 and A 2 7^ 1 has the phase diagram qualitatively similar to that shown in Fig.l [24| . 
The point |J 2 /Ji| = 1/4 is the transition point for the isotropic model (Ai = A 2 = 1), where the transition from the 
ferromagnetic to the incommensurate ground state with S z = occurs. This transition is the second order one. For 
the anisotropic model the phase transition between the ferromagnetic and the incommensurate states is of the first 
order type. 

The properties of the commensurate spin- fluid and the incommensurate phases have been studied in [24]. In 
this paper we study the properties of the ferromagnetic (F) phase. Though the ferromagnetic ground state is very 
simple, the excitation spectrum is not trivial. The important feature of the spectrum is the existence of two types 
of excitations: conventional spin waves and the multi-magnon bound states. These excitations govern the low- 
temperature thermodynamics of the model in the F phase. The additional motivation of the study of this phase is the 
fact that the ferromagnetic in-chain ordering has been observed in some edge-sharing cuprates [25| . This indicates 
that the frustration parameter A = | J 2 /Ji| can be less than 1/4 in these compounds [261 ]. 

In our study we focus on the excitation spectrum and the low-temperature thermodynamics of the F-AF anisotropic 
model in the F phase. For simplicity we confine our study to the case of the anisotropy of the NN interaction only, 
i.e. we consider Hamiltonian ([1]) with Ai = A > 1 and A 2 = 1 

H = - y^(5^5^ +1 + Sn&n+l + ASn&n+l ~ 4") + ^^(^n " &n+2 - ^) (2) 
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FIG. 1: The phase diagram of model © 



where we take |Ji| as an energy unit and add constant shifts to secure the energy of the fully polarized state to be 
zero. 

We are mainly interested in the behavior of the model in the vicinity of the transition point (A = 1/4, A = 1). 

The paper is organized as follows. In Sec. II we represent the known results for model ^ at A = 0. In Sec. Ill 
we exactly calculate the two-magnon bound energy. In Sec. IV we perform the scaling estimates for the energy of 
the multi-magnon bound states based on the analysis of infrared divergencies in the perturbation theory in small 
parameter a = A — 1. In Sec.V we present results of numerical calculations of finite chains. In Sec. VI we study 
the low-temperature thermodynamics of model ^ and determine the region of parameters where the multi-magnon 
excitations dominate. The relevance of the F-AF anisotropic model to the copper oxide L%iCuO/± is discussed in 
Sec.VII. 



II. BOUND STATES IN THE FERROMAGNETIC CHAIN 



In this section we review the known results relevant to our study. In the special case A = model ((2]) reduces 
to the ferromagnetic XXZ chain with the Ising-like anisotropy. Its ground state is ferromagnetic with zero energy 
and with the gap in the excitation spectrum. Lowest-lying excitations are bound states of overturned spins from the 
fully polarized ground state (multi-magnon bound states). The energy of m-magnon bound state E m (k) for the chain 
with the periodic boundary conditions (ring) was found by Ovchinnikov (2?| using the Bethe Ansatz. He showed that 
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E m (k) at N — > oo is 



,, , sinh (z/) [cosh(mzy) — cosfcl 

£™ K = r-TT ^ (3) 



where A: is a total momentum, cosh^ = A and m = 1,2, . 
For fc = the above equation reduces to 



£ m = E m (0) = sinh (i/) tanh (^J (4) 



In particular, for one- and two-magnon states it gives 

E 1 = A - 1 

E 2 = A-~ (5) 

(certainly, there is no bound state for m = 1, in this case Eq.([3]) describes the one-magnon spectrum). 

It follows from Eq.([3]) that the energy of the m-magnon bound state saturates exponentially with m to the value 
E s = VA 2 — 1 and does not depend on k for m>l, i.e. excitations become dispersionless. 

Another important result for model ([2]) at A = was obtained by Alcaraz et.al. in Ref.[28|, where the multi-magnon 
bound states were studied for the chain with free-end boundary conditions (open chains) . The Bethe Ansatz solution 
gives the energy of the m-magnon bound state in the open chain in a form [281 ] 

m 1 1 

E m , open = ^2(A--( Pi --)) (6) 



where pi is defined by a recurrence equation 



2A - p t 



with pi — 1/A. 

It can be shown [2^] that the solution of Eq.([7|) gives for -E mj0 pen 

-Em.opcn = i^A 2 - ltanh (mv) (8) 
The comparison of Eq. ([8]) with Eq. ([4} leads to the exact relation 

^2m — 2_£/ m0 p Cn (9) 

This relation means that the energy of the multi-magnon bound state in the open chain at m 3> 1 is a half of that 
for the ring. This property has been observed earlier [30( in numerical calculations of model @ for A = 0. Actually, 
the validation of Eq. ([9]) looks quite natural, because it implies that the magnetic soliton of size 2m can be represented 
as two kinks of size m. 

As was shown in (30j a low-temperature thermodynamics of the anisotropic ferromagnetic chain is determined by 
an effective gap which is the lowest of two values: the gap for the one-magnon excitations E\ (spin-waves) and the 
gap for the multi-magnon bound states for the open chain (not ring) E m ^ opcn . The comparison of Eq.([5]) and Eq.([8]) 
shows that the effective gap for A < 5/3 is the spin wave gap E\ = A — 1, while for A > 5/3 it is the multi-magnon 
bound state energy equal to -Ekink = |VA 2 — 1. 

As will be shown below, many peculiarities of the excitation spectrum of the anisotropic ferromagnet remain for 
model © with A ^ 0. 

III. ONE- AND TWO-MAGNON EXCITATIONS 

In this Section we consider the one- and two-magnon states for the frustrated model @ with A ^ 0. We begin with 
the case of the rings with periodic boundary conditions. The energy of the one-magnon state is 



Ei(k) = A-cosfc- A(l-cos(2fc)) 



(10) 
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The one-magnon spectrum has a minimum at fc = for < A < 1/4 and has a double- well form with two minima 
at k = ± arccos (1/(4A)) for A > 1/4. The expansion of Ei(k) at small k (a = A — 1) 

£i(fc) = «+ 1 ^ 2 + ^ r ^ 4 + ... (ID 

shows that the behavior of the low- lying excitations on the isotropic line A = 1 is different: in the region < A < 1/4 
the low-lying excitations are described by k 2 spectrum 

E 1 (k)= l —^k\ A<l/4 (12) 
while at the transition point A = 1/4 the one-magnon spectrum becomes 

E 1 (k) = !j, A =1/4 (13) 

As was shown in Ref. [24( this difference plays a key role in the change of the critical exponents near the transition 
point A = 1/4. 

A remarkable feature of the two-magnon spectrum is the existence of the bound states lying below the scattering 
continuum. For the isotropic F-AF model (A = 1) the two-magnon bound states at A > 1/4 have been studied in 
detail earlier 

i, US [HI, Hi!, M, i| . In particular, it was shown that namely these states define the saturation magnetic 
field in the incommensurate phase for A > 0.367 [34| . 

Fortunately, the two-magnon bound state energy E2(k) can be found exactly for general anisotropic case of model 
([2]) . The analysis of the scattering problem of two magnons shows that for each total momentum k of the magnon pair 
there is one bound state. The minimization of the energy E2{k) over k gives the gap in the two-magnon spectrum. 
In general, the dependence of /c m ; n (A, A) minimizing E2(k) is rather complicated. We note only that fc m i n — * tt when 
A is increased at fixed A and the dependence of the value A^ (A) was found in Ref. [3l[ . 

We are interested mainly in two-magnon bound states in the ferromagnetic phase in the vicinity of the transition 
point A = 1/4 and for weak anisotropy a<l, In this region the energy E2{k) has a minimum at k = 0. Therefore, 
hereinafter we restrict ourselves to the magnon excitations with k — 0. For this case the two-magnon bound state 
energy E2 is determined from the equation: 

1 } E 1 (q)dq , . 

1 1 ' = -1 (14) 



tt J E 2 -2E 1 (q) 


where E\ (q) is the spectrum of the one-magnon excitations given by Eq. (jTT 

Evaluating the integral in Eq. (fT4"|) we obtain the following algebraic equation for E2 



[(2A~£ 2 ) 2 -4] 1 - 2A(2A - E 2 ) + 2A^(2A - E 2 ) 2 - 4 = E 2 (15) 

For A = the solution of this equation reproduces the exact result Eq.(j5|). 
Let us represent E2 in the form 

E 2 = 2a-^Q{- 1 )-E b (16) 

where is the Heaviside function, 7 = 1 — 4A and E b is the binding energy relative to the two-magnon scattering 
continuum. 

The series expansion of the binding energy in small parameter a for A < 1/4 is: 

E b = — -4H^ + l 1/2 -l) + 0(a i ) (17) 
At the transition point A = 1/4 the solution of Eq. (fT5l) gives 

E b = a 4 / 3 - ^a 5 /3 + (a 2 ) (18) 
Near the transition point (a <C 1, |-y | <C 1) Eq. (fT5)) can be simplified to the cubic equation 

(E b + 7 2 e(- 7 )) 3/2 + i{E b + 7 2 e(- 7 )) = « 2 (19) 
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The solution of Eq. (fT9| has a scaling form depending on a scaling variable k — 7/a 2 / 3 

E b = a^ 3 g( K ) (20) 

where the function g(n) is a solution of equation 

(g + k 2 0(-k)) 3/2 + K (g + k 2 Q(-k)) = 1 (21) 

In the limits k — > 00 (7 3> a 2 / 3 ) the asymptotic of g(n) reproduces the leading terms of Eq. (fT7|) . In the region close 
to the line A = 1/4, when j <§; a 2 / 3 (\k\ <C 1) the expansion of g(n) = 1 — 2k/3 + . . . results in the following series 
for the binding energy 

E b = a 473 - ^a 2 ? 3 + 0( 7 2 ) (22) 



which contains the corrections to the leading term in Eq.fTSj). 

In the region A ^ \ the binding energy vanishes on the curve 2a — j 2 . Near this curve according to Eq. (|15[) the 
binding energy behaves as 

Eb - 2| 7 | (23) 

Thus, the critical exponent characterizing a power law dependence of the binding energy on a changes from 2 at 
A = to 4/3 at A = 1/4. This change is due to the modification of the behavior of the one-magnon energy at small k 

Eqs. cental. 

For further study it is useful to estimate the finite-size corrections to the two-magnon binding energy. For finite 
rings the integral in equation (114[) is replaced by a sum: 

l y _#L = 

N^E 2 - 2E 1 (k) ~ [M) 

k 

with jfc = 2™ and n = -N/2, . . . N/2. 

The leading terms for the binding energy are governed by small k, when denominators in Eq.(|24p are small. 
Therefore, we take into account only the leading term in the expansion of E±(k) fEq. (|lip ) and take infinite limits of 
sum in Eq.([2I|). So, for A < 1/4 Eg. (|2"3|) reduces to 

2a ^ 1 

— > 2 = 1 25 



n— — 00 &b ^~ 7 V TV 

Evaluating the sum we obtain the following equation for 



a cothf^)=l (26) 



VWb V 2^7 

At N ^> j/a one obtains the exponentially small finite-size correction for the leading term of the binding energy 

2 

E b = — (\ + Ae~ Nah ) (27) 
7 

Similar procedure for the case A = 1/4 results in the equation 

2a °° 1 

~~N £ Sh+ W2^ ) 4= 1 ( 28 ) 

which leads to the following equation for E b 

a sinh (iV^ 174 ) + sin (nE^) 
_ ^ L ^ L = 1 (291 

E 3 b /4 cosh (TVS 174 )- cos (NE^) 
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At Not 1 ' 3 > 1 it gives for the binding energy 




Thus, we found the exponentially small finite-size effects, which is not a surprise for the bound states. More 
important is that we identified the scaling parameters iVa/7 for A < 1/4 and Na 1 / 3 for A = 1/4, which will be 
exploited later. 

At the end of this section it is worth to make a following remark. Let us consider the one-magnon states in the open 
chains. The one-magnon problem in this case can be solved by a standard method. The spectrum of the one-magnon 
excitations consists of (N — 2) band states, the energies of which coincide at N — > 00 with those given by Eq. pi)]) . 
However, there are two degenerated bound states localized near chain ends. Remarkably, the energy Ex t0pen of these 
bound states is determined by Eq. (|15[) with the substitution of 2_Ei i0pcn for E<x, Thus, the energy of the one-magnon 
bound states in the open chain is half of the energy of the two-magnon bound state with k = in the ring. It means 
that relation @ remains valid for the case m = 1. 



IV. PERTURBATION THEORY 



For more than two magnons the exact analytic solution is not possible excluding the case A = 0. Therefore, for 
studying the multi-magnon bound states we develop the perturbation theory (PT) in small parameter a. 

At first, we inspect how the energy of the two-magnon state obtained in the preceding section can be estimated in 
the framework of the PT. 

A. Perturbation theory for two-magnon states 

Let us represent Hamiltonian (|2jl in a form 

H = H 1 +XV 2 +aV z 
Hi 



V 2 



- ^(Sn • S„+i - -) 
y^(Sn • S„ +2 - -) 

-T,( S nS Z n+l-\) (31) 



We use the perturbation theory in small parameter a, so that a small perturbation aV z is added to the unperturbed 
Hamiltonian Hq = Hi + AV2. The ground state of Ho for A < 1/4 is ferromagnetic (total spin S — N/2) and is 
degenerate with respect to total S z . The perturbation aV z splits this degeneracy. We consider the two-magnon 
sector, therefore, we develop the PT to the ferromagnetic state with the projection S z = N/2 — 2: 

M =7m^T) {S ~ rm --- n (32> 

where S~ = ^2 S~ is the total descending operator. 
In the first order of PT we have 

fi (1) = (MaV z \%) =2a (33) 

The second order reads 

E {2) = a 2 (^o\Vz\^ q ) (jj q \Vz\ipo} , 34 s 

The state |^o) has the total momentum k = 0. The perturbation V z preserves the total momentum k and the 
projection S z , but can change the total spin S on value ±2. Therefore, only two-magnon states \ip q ) with k = and 
S = N/2 - 2 are involved into the PT 

i^)=E^( / )^s- +i irr...t) (35) 

nj 



7 



where norm (ip q \i> q ) = 1 requires (f q (l) ~ 1/iV. 
The matrix elements are 



(36) 



The iV-behavior of denominators in Eq. (|34p depends on the value of A. For A < 1/4, the one-magnon spectrum is 
Ek — r yk 2 /2 fEq. ([T2")) ) and denominator behaves as 

E q -E ~ 1 N- 2 (37) 

For A = 1/4, the spectrum is e k = k 4 /8 fEa.([T3))) and 

E q -E ~ N- 4 (38) 

Thus, the estimate of the second order Eq. (|3~4"|) gives 



E^ ~ - — , A<l/4 
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a 2 N 2 , A = 1/4 (39) 

From Eq. ([39|) we see that the second-order correction for the ground-state energy converges for A < 1/4, but 
diverges for A = 1/4. Therefore, we need to estimate the higher orders of the PT. The perturbation series for the 
two-magnon binding energy fEq. (|16p ) can be written in a form: 

E b (a) = - (Vol a 2 V z 1 V z + a 3 V z 1 V z 1 V z + ... \^ ) (40) 

&n — tin rjo, — tlo, tin — Hn 



Suppose that the main contributions to the energy are given by the low- lying excitations Eqs. ((37|l . (|38|) . Then the 
higher orders of the perturbation series contain more dangerous denominators and, therefore, possibly have higher 
powers of the infrared divergency. Therefore, we use scaling arguments to estimate the critical exponent for the 
ground-state energy. Below we will pay attention only to the powers of the divergent terms and omit numerical 
factors. 

The analysis shows that the matrix elements of the perturbation operator V z between the states \ip q ) involved in 
the PT behave as 

m V z |Vy) = AW1W(1) ~ l/# (41) 

Collecting the most divergent parts in all orders of the PT, we express the correction to the binding energy as: 

E b = a (ip g \ V z |Vy) Y, c nx n = ^f(x) (42) 

where c n are unknown constants and 

x - ^A^kl (43 ) 

Eg — Eq 

is a scaling parameter, which absorbs the infrared divergencies. 
The estimate of scaling parameter using Eqs. ([3^)) - ((551) gives 

x=^, A<l/4 
7 

x = aN 3 , A = 1/4 (44) 

We note, that exactly the same scaling parameters were determined in previous Section fEqs. ([27|) . (|30p ). 
In the thermodynamic limit the binding energy tends to a finite value (independent of AT). Therefore, the scaling 
function at N — * oo behaves as 

/Or) ~ x, A<l/4 

f(x) ~ x 1 / 3 , A = 1/4 (45) 
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This gives for E b : 

a 2 

E b , A < 1/4 

7 

E b ~ a 4 / 3 , A = 1/4 (46) 

The obtained results totally agree with those obtained in previous Section. Thus, the scaling estimates of diver- 
gencies in the PT correctly reproduces the scaling parameters and the leading terms in the energy. 

B. Perturbation theory for m-magnon states 

1. Perturbation theory at A = 

We start to study the multi-magnon problem from the exactly solvable case A = 0. The lowest m-magnon energies 
in the thermodynamic limit are given by Eq.(j4])- For small anisotropy a <C 1 the energy E m can be written as 

E m {a) = maf (am 2 ) (47) 

with 

f (x) = \/^tanh ( x ff) (48) 



x VV 2 , 

So, the expansion of Eq. (|47[) in small parameter am 2 is 

E m (a) = ma - -m 3 a 2 + —m 5 a 3 + . . . (49) 
v ' 6 30 v ; 

However, in the opposed limit am 2 3> 1 related to large magnon complexes, the energy converges to 

E s {a) = \/2a (50) 

On the other hand, from the solution of two-magnon problem we know that for m = 2 but hnite N , the energy 
depends on the scaling parameter x — aN (Eq. (f27|) or Eq. ((42|) ). However, for m > 2 this scaling parameter is modified 
and becomes x — amN 29!]. Thus, for the case when both m and N are large but finite the energy has a scaling form 
of two scaling parameters: 

E m (a) = mafi(am 2 ,amN) (51) 

with fi(am 2 , amN) — > fo(am 2 ) at N — > oo. 



2. Perturbation theory for < A < 1/4 

Now let us consider the case A -C 1. In this case PT pip contains two small parameters a and A and, consequently, 
two channels V z and V2 ■ Each channel can produce infrared divergencies and is described by its own scaling parameter 
[35] ] . The matrix elements of the operator V2 at A <C 1 were found in Ref. [lij : 

(Vi|Wi)~W- 2 (52) 

Such behavior of the matrix elements means that the perturbation V2 does not produce infrared divergencies and, 
therefore, does not form the scaling parameter [24|. It is natural to expect that the behavior of the matrix elements 
of type (|52p remains the same up to the point A = 1/4. This assumption results in the following expression for the 
lowest energy in the sectors with small value of S z (am 2 ^> 1) in the region < A < 1/4: 

E s (a,X) = A(\)Va~ (53) 

where A(X) is a smooth function with A(0) = a/2 (see Eq.(j5U|). 
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3. Perturbation theory at A = 1/4 



At the transition point A = 1/4 the Hamiltonian has a form 



H = Hx + -V 2 + aV z 



(54) 



In order to make the scaling estimates of the PT in a one needs to know the matrix elements of the operator V z . 
Unfortunately, it is very difficult problem. However, we can assume that the behavior of matrix elements with m 
and A for A = 1/4 remains the same as for the case A = 0. Then, the only difference between these two cases is the 
modification of the spectrum presenting in denominator in Eq. |43|) . According to Eqs. (|37|) .([38 |) for the case A = 1/4 
both scaling parameters in Eq. ([5l"j) acquire additional factor A 2 . So, for A = 1/4 Eq. ([5T]) transforms to 



E m {a) = maf xy (x,y) 



(55) 



with x = amN 3 and y = am 2 N 2 . 

In the thermodynamic limit, when both x — > oo and y — > oo, the scaling function f X y(x,y) becomes a function of 
one variable (independent of A) 



y=^ = a 1/3 m 4 / 3 

x I 

and the m-magnon energy takes the form 

E m (a) = maf v (v) 
We expect that the expansion of the scaling function f v {v) at v = is 

f v (v) = 1-C x v + C 2 v 2 + ... 
which gives the following expansion for the energy: 

E m (a) =ma- Cim 7/3 a 4/3 + C 2 m 11/3 a 5/3 + 



(56) 

(57) 
(58) 
(59) 



So for the case to = 2 it reduces to the obtained result for the two-magnon binding energy Eq.fTSj). 

For sectors with small value of S z (m ~ A/2), the scaling parameter v ^> 1 and the energy E m (a) saturates to the 
finite value E s (a) describing the finite gap in the spectrum. Therefore, the asymptotic of the scaling function j v (y) 
at v — ► oo is fviy) ~ v~ 3 l^ and, finally 



E s (a) = Ba 3/4 



(60) 



where a constant B will be found numerically. 



4- Perturbation theory near the transition point 
In the vicinity of the transition point the PT contains two perturbations aV z and JV2 

H = H 1 + -V 2 + aV z - ^-V 2 (61) 
4 4 

We assume that the behavior of the matrix elements of the operator V 2 remains the same as in the region A < 1/4 
(see Eq. (j52"10 . At A = 1/4 according to Eq. (|4*5)l the perturbation ^V 2 produces the scaling parameter z — 7 A 2 24]. 
Thus, near the transition point A = 1/4 the TO-magnon energy can be written in a scaling form 

E m (a, 7) = maf xyz {x, y, z) (62) 

In the thermodynamic limit, when all parameters x,y, z — > 00, the scaling function f X yz{x, y, z) becomes a function 
of two independent of A variables: v fEq . ([SB"]) ) and 

/i = = —= 63) 

x Ja 
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It is worth to note that the scaling parameter k obtained in the exact solution of two-magnon problem near the 
transition point Eq. (|20p transforms to the parameter /i Eq. (|63[) with increase of m. The m-magnon energy in the 
thermodynamic limit takes a form 

E m (a, 7) = ma/^(/i, u) (64) 

For the lowest states in sectors with small value of S" (m ~ N/2), the scaling parameter v 1 and the multi-magnon 
energy becomes 

E s {a n ) = o?l i f fl {n) (65) 

The function / M (/Lt) is generally unknown, however we can determine its behavior in two limits. Close to the line 
A = 1/4, when \i — > 0, Eq. (p5|) reduces to Eq. (|6"0)) . which implies that /^(O) = B. In the limit /i — > 00 Eq. ([631> must 
reduce to Eq. (|53[) . This requires the asymptotic /^(/x) ~ Vm at — > 00 and the energy 

£■3(0,7)^^/70, a < 7 2 (66) 

This in turn means that the function A(X) in Eq. ([53|) behaves at A — + 1/4 as 

A(X) ~ VI -4A (67) 



V. NUMERICAL CALCULATIONS 



We have carried out exact diagonalizations of finite rings and open chains up to 24 sites. We observed that the 
multi-magnon bound states are formed when the corresponding finite-size scaling parameters becomes large. This 
imposed the natural restrictions on our calculations: a 3> 1/jV 2 for A < 1/4 and a ^> 1/N A near the transition 
point. We found that in all considered cases the multi-magnon bound energies have exponentially small finite-size 
corrections (for two-magnon binding energy this fact was established analytically in Sec. II). Therefore, we used a 
linear extrapolation in exp(— aN) with a fitting parameter a. To check the numerical accuracy of the extrapolation 
we compared the extrapolated energy for the total S z = (m = N/2) and A = with Eqs.(|4]) and ||5J) and found 
perfect consistency with the exact results. Thus, for not too small a (a ^S> 1/N A near the transition point) the use 
of finite chains with N < 24 is sufficient for the extrapolation to the thermodynamic limit. We note also that the 
saturation of the energy E m at m — > N/2 and the convergence to the thermodynamic limit for the open chains occur 
noticeably faster in comparison with the rings. Really, the kink excitation on open chain of length N corresponds to 
the soliton excitation on ring of length 2N . Therefore, in our numerical calculations we used mostly the open chains 
in subspace S z = 0. 

Firstly, we validated the important relation ©. In Section II we have shown that for A = and k = the energy 
of the 2m -magnon bound state on the ring is double of that for the m-magnon bound state on the open chain. We 
showed also that at N — > 00 relation @ is valid for m — 1 and A ^ 0. Unfortunately, we can not rigorously prove 
this relation in the general case when m > 1 and A / 0, because the analytic solution is not possible in this case. 
However, we checked numerically that the relative difference (2i? mj0 p Cn — E2 m ) / ' E% m vanishes rapidly with N, so that 
for N = 24 this difference is less than 0.1% for m = 1 ... 6 in a wide range of values of a and A < 1/4. Therefore, 
we suggest that at N — > 00 the energies 2£ , m , opcn and i?2m do coincide. The typical dependence of the bound state 
energy on m for rings and open chains of length N — 24 is demonstrated in Fig. 2. As follows from Fig. 2 both energies 
2£ m ,opcn and E m saturate to the same finite value when m 3> 1. Hence, we expect that in the thermodynamic limit 
and for m ~ N/2 the energy of the magnetic soliton is double the kink energy 

E s (a,X) =2E kiDk (a,\) (68) 

In Fig. 3 we demonstrate the dependence of the extrapolated lowest energy -Ekink(a) in the sector with S z = for 
A = 1/4. A linear fit of E k i nk /a 3 ^ 4 as a function of a 1 / 4 gives the equation 

E kink ~ 0.35a 3/4 + 0.275a (69) 

As one can see the correction term 0.275a gives substantial contribution for not too small a and can not be neglected 
(for example for a = 0.1 it gives near 30% of the kink energy). Therefore, one has to take it into account and rewrite 
Eq.® as 



(70) 
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m 



FIG. 2: Energies of the m-magnon bound states of model ([2} for rings (circles) and open chains (triangles). Arrows connect a 
few pairs of points corresponding to double m-magnon energy for open chain 2_E , mi o pC n and 2m-magnon energy for ring Ei m , 
validating Eq.(|9}. 



The linear terms in the expansion of functions ffj,(n) and in /j, correspond to the first order of PT (|61|) in 

%V-2 and are defined by the behavior of the correlator (S„ • S n+ 2 — j) at A = 1/4. The extrapolation of numerical 
calculations gives 

]T(S„ • S„ +2 - ~) ~ 1.92a 1 / 4 - 0.84a 1 / 2 (71) 

n 

This means that for \fx\ <C 1 the functions and are 

/„(/*) ~ 0.35 + 0.48^ 

g^{p) ~ 0.275- 0.21^ (72) 

It turns out that near the transition point it is sufficient to take the function in a form of Eq. ([72^) . This 

fact is confirmed by Fig. 4, where we present the dependence of fn(^) — (-Ekink — Q 3At(/ i ))/ Q;3 ^ 4 as a function of the 
scaling variable [i = j/y/a, calculated for different values of a and 7 in the range a — ... 1 and 7 = —0.2 . . .0.2. 
All calculated points lie perfectly on one curve (see Fig. 4). This confirms the scaling form (fT0"|) in the vicinity of the 
transition point. As we see in Fig. 4 at <C 1 the function has the expansion (|72|) . In the limit n ^> 1 the 

numerical calculations gives /^(/i) — 0.65^/jEt. 

To end this section we list the main results obtained by numerical calculations of finite chains for A ^ 0. The 
m-magnon bound state energy saturates for m 3> 1 for both rings and open chains, describing the finite gap in 
the spectrum. The multi-magnon bound complexes with m 3> 1 are very massive, which results in flat band and 
quasi-degeneracy of the bound m-magnon excitations for the rings over the total momentum k. This means that the 
total degeneracy of the 'soliton' energy level E s for the rings is proportional to N 2 . In contrast to the rings, the 
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FIG. 3: The dependence S killk /a 3/4 vs. a 1/4 for open chains at A = 1/4. Linear fit corresponds to -Bkink = 0.35a 3 / 4 + 0.275a. 

m-magnon bound state energy for the open chains is two-fold degenerated for each m. Therefore, for the open chains 
the degeneracy of the kink energy level -Ekink is proportional to N. The kink and soliton energies satisfy relation 
(I68[) . These properties resemble those found for A = 0, though the energy of multi-magnon excitations is modified for 
A 7^ 0. In particular, the critical exponent depends on A. 

It may be of interest to remark one point related to the bound magnon states in open chains. The lowest m-magnon 
state is the state with one domain wall (kink state) and it describes the gapped excitation above the fully polarized 
ground state. However, it is possible to realize this kink as the ground state. Let us add to Hamiltonian @ the term 
with boundary magnetic fields 

- HS{ - (73) 

Then there is a critical magnetic field h c for which the kink energy is zero and the kink state is degenerate with 
the ferromagnetic state. It was shown in |28| that h c = i\/ A 2 — 1 for A = and the critical field h c does not depend 
on m. At h = h c the ground state is iV-fold degenerate. It is a result of the special symmetry of Hamiltonian ((2J) 
for A = with the boundary magnetic field h = h c . This is not the case for A ^ 0. In this case the magnetic field 
h c depends on m. For example, h c (l) = a 3 / 4 /2 1 / 4 ~ 0.84a 3 / 4 for a <C 1 and A = 1/4. The value h c (m) decreases 
with m and saturates to h c ~ 0.35a 3 / 4 for m 3> 1 and A = 1/4. For h > h c the kink state with m — N/2 (S z = 0) 
becomes the ground state. 

VI. LOW-TEMPERATURE THERMODYNAMICS 

Results of numerical calculations show that many peculiarities of the low-energy spectrum of model © at A ^ 
are similar to those for the anisotropic ferromagnetic chain. The thermodynamics of the latter model was studied in 
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FIG. 4: The scaling function ffi(fi) in Eq. (|65[l for open chains. Circles correspond to different small values of 7 and a. Circles 
form perfectly one curve, justifying the scaling dependence (|65ll . 



Ref. [3C| and we use the arguments of this work to treat the low-temperature thermodynamics of model © . As was 
shown in Ref. [30j| the principal contributions to the partition function at low temperatures are given by two classes 
of the excited states: spin waves and the multi-magnon bound states. Both types of excitations are gapped and the 
existence of the gap implies an exponential behavior of the thermodynamic functions at T — ► 0. 

The leading terms for the free-energy at T — ► are given by a sum over low- lying states of both types. Certainly, 
the free energy must satisfy the natural condition to be proportional to N. This requires that the number of the 
excited states of each type must be proportional to N too. As for the spin waves the number of such states is ~ N for 
both rings and open chains. But it is not the case for the multi-magnon bound states. As was shown in the preceding 
section the number of these states for open chains is ~ N, while for rings it is ~ N 2 . Therefore, as was already 
noticed in Sec. II, the correct contribution of the multi-magnon excitations to the free energy is given by those for the 
open chains (kinks) rather than for the rings (solitons): 



Fbu 



-NTexp 



Ekinkja, A) 
T 



(74) 



where we put Boltzmann's constant ks = 1. 

Of course, there is no difference between open chains and rings in the thermodynamic limit and Eq. (|74p can be 
obtained using the rings as well. In this case a more complicated summation of high-lying excitations leads to the 
effective double reduction of the soliton energy E s , which restores Eg. ([74")) . 

The spin-wave contribution F sw to the free energy is 



F — 



NT 



exp 



E^k) 
T 



dk 



(75) 
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where Ex(k) is one-magnon energy (|10p . 

Evaluating the integral in Eq. ([75]) we obtain for A < 1/4 



F sw = -^^- a) ' T K l/4 ( 7 2 /T) (76) 



where K\u is the modified Bessel function of the second kind. 
In particular, for T <C J 2 



For A = 1/4 (7 = 0) 



For A > 1/4 and T< 7 2 



NT 3/2 - a /T 

(77) 



/2-7T7 



Fsw = ~ 23/4r(3/4) (78) 



NT 3/2 y 2 /2T-a/T 

F sw = - ' (79) 

V 71 " ItI 

The dominant contribution to the low-temperature free-energy is given by the excited states with the minimal value 
of the gap. Therefore, in order to identify the prevailing type of excitations we need to compare the kink energy -Ekink 
with the spin wave gap E sw , which corresponds to the minimum of one-magnon spectrum (|10p 

E sw = a, A < - 

^sw = «~ y. ( 8 °) 

It turns out that there are two regions in the phase plane (A, a) where the low-temperature thermodynamics is 
governed by different excitations (see Fig. 5). The boundary between these two regions a c (A) is determined by the 
equation 

-Bkink(a c , A) = E sw (a c , A) (81) 

The calculated dependence a c (A) is shown in Fig. 5. In the region of the F phase below the curve a c (A) the dominant 
excitations are spin-waves, while for a > a c (A) they are multi-magnon bound states. As follows from Fig. 5 the value 
a c (A) decreases when A increases. In particular, a c (0) = 2/3 and a c (l/4) ~ 0.046. 

The specific heat in the region with the spin-wave dominance is 

c = ^hm' T<<7 

a 2 e- a l T 1 
T 7/4 2 3/4 r ( 3 /4)' ~ 4 

C = [ 21 - , A>i T« 7 2 (82) 

t 3/2 7^M ~ 4 

In the region a > a c where the dominant contribution to the low-temperature thermodynamics is given by the 
multi-magnon bound excitations the specific heat has an Ising-like behavior 

C=^exp(-^) (83) 

To obtain the low-temperature susceptibility we add the magnetic field h along the Z axis. Then the spin-wave 
free energy F sw has the same form as given by Eqs. ([75|) with a replaced by (a + h) and the spin- wave contribution to 
the susceptibility is 

X S w = -^f (84) 
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FIG. 5: The dependence a c (X). For a < a c 
a > a c by multi-magnon bound states. 



the low-temperature thermodynamics is governed by spin-wave excitations, for 



In the presence of magnetic field the energy of m-magnon excitations is E rn + mh. The contribution to the low- 
temperature susceptibility of the multi-magnon bound excitations can be obtained using the following arguments 
[3(il ] . At ra > 1 the magnons are tightly bound in the multi-magnon complex and the size of the m-magnon bound 
state is m with the exponentially accuracy in m. Besides, these complexes are immobile because of their very 
large mass. Therefore, these bound states can be considered as domains of m neighbor overturned spins in the 
one-dimensional Ising model with the effective exchange constant £ ; kink(a, A). Using this analogy we obtain the 
multi-magnon contribution to the zero-field susceptibility in the form 



1 ( -Ekink 

exp 



XT 



(85) 



Comparing Eq. (|84p and Eq. (|85p one can easily see that the multi-magnon contribution to the low-temperature 
susceptibility is dominant at all a and A in contrast to the case of the specific heat. 

We note that equation (|85|) determines the low-temperature susceptibility for the anisotropic model, i.e. for A > 1 
and it is not valid in the limit A = 1, because the symmetry of the Hamiltonian changes at A = 1. For example, the 
susceptibility of the isotropic ferromagnetic Heisenberg chain (A = 1, A = 0) x ~ 1/T 2 . The susceptibility for the 
isotropic F-AF model can be found using the modified spin wave approach (MSWT) [36]. In particular, for A = 1/4 
and A = 1 this method gives [2!| 



1 



Xniswt 



4(2T) 4 /3 



(86) 
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VII. DISCUSSION AND SUMMARY 

Now, let us discuss a relevance of the considered model to the copper oxide Li 2 CuO^. This cuprate consists of 
the chains formed by the edge-sharing CuO^ squares (37[- The magnetic interaction between NN spin-1/2 Cu 2+ 
ions along the chain is ferromagnetic while an exchange interaction between NNN Cu ions is antiferromagnetic and 
these chains are described by the F-AF model. The magnetic structure of Li 2 CuO^ was determined by neutron- 
scattering experiments (25[. Below the Neel temperature ~ 9K the spins of each Cu0 2 chain has a ferromagnetic 
arrangement and the arrangement between neighboring chains is antiferromagnetic. At present the reason for the 
observed ferromagnetic in-chain order is unclear. The early estimations of the frustration parameter gave A ~ 0.4 — 0.6 
[lH [37j . But for such values of the frustration parameter model ([2]) has a spiral-like ground state rather than the 
ferromagnetic one [5j, l8| . To resolve this discrepancy it was proposed that the observed ferromagnetic order arises 
due to the specific role of the interchain interactions [53, [Hj]. However, recent estimations [26| of this value based on 
both exact diagonalization of Cu — O Hubbard model and the DFT calculations show that the frustration parameter 
is somewhat smaller than the critical value (A ~ 0.23). Therefore, we can suppose that individual, non-interacting 
chains have the ferromagnetic ground state. The AF long range order in Li 2 CuO± below the Neel temperature arises 
due to a weak antiferromagnetic interchain interaction J± . 

A standard method for treating the quasi-one-dimensional systems is the mean-field approximation for the interchain 
interaction. In this approximation the Neel temperature is determined by the equation 

zJ± X id(T n ) = 1 (87) 

where xid is the susceptibility of the individual chain and z is the transverse coordination number. In Li 2 CuOi each 
Cu0 2 chain is surrounded by four parallel neighboring chains and the only non-negligible coupling J± occurs between 
NNN spins on neighboring chains. Therefore, the effective coordination number is z = 8. 

To estimate T N from Eq. ([57|) we use the susceptibility Xid(T) given by Eq. (f8"5]) . According to Ref.[26| the frustration 
parameter for Li 2 CuO^ is close to the critical value and we take for the gap -Ekink^, A) its value at A = 1/4 and a« 1: 
£kink = I Jil(0.35a 3 / 4 + 0.275a). The NN in-chain interaction J\ and the interchain interaction J± were estimated in 
Refs.(26l. [39| as J\ ~ — 145^ and J± ~ 3.6K. Using these parameters we calculated the dependence Tjv(a) shown 
in Fig. 6. According to these calculations the anisotropy in this compound is estimated as A ~ 1.01. Certainly, 
this estimate is based on the mean-field treatment which overestimates the transition temperature. Nevertheless, 
we expect that the anisotropy of the exchange interactions in Li 2 CuO^ does not exceed a few percent. This fact is 
also confirmed by the estimate of the Neel temperature for the isotropic case (see Fig. 6), where we used the MSWT 
susceptibility (|56")1 in Eq. ([57|) . The so obtained Neel temperature Tjy — 7. 5 A" is slightly lower than the experimental 
value Tjv — 9A, confirming the presence of weak anisotropy. Though the anisotropy is very small it can essentially 
affect the excitation spectrum especially when the frustration parameter is close to the critical value A = 1/4. 

In conclusion, we study the excitation spectrum of the one-dimensional anisotropic F-AF model in a parameter 
range corresponding to the ferromagnetic ground state. The remarkable feature of the spectrum is the existence of 
the multi-magnon bound states. The lowest-lying m-magnon excitations are quasi-degenerated and are separated by 
the gap from the ferromagnetic ground state. This gap as a function of a at a <C 1 has a power-law behavior with the 
exponent depending on the frustration parameter. It turns out that the gap for the bound multi-magnon excitations 
in the rings is twice as large of that in the open chain. The multi-magnon excitations together with the spin waves 
give dominant contributions to the low-temperature specific heat. The thermal gap characterizing the exponential 
behavior of the thermodynamic functions is the smallest value of two gaps: for one-magnon excitations and for bound 
multi-magnon ones in open chains. The comparison of these gaps defines the regions of the dominance of one or 
another type of excitations. Contrary to the specific heat the zero-field susceptibility is always determined by the 
multi-magnon excitations. 
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